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Abstract 

We discuss cosmology of models with universal extra dimensions, where the Standard Model 
degrees of freedom live in a 4 + n dimensional brane, with n compact and small extra spatial 
dimensions. In these models, the simplest way to obtain the conventional 4-dimensional Planck 
scale starting with a low string scale is to have also some larger extra dimensions, where only 
gravity propagates. In such theories, dimensional reduction generically leads to at least two radion 
fields, one associated with the total volume of the extra spatial dimensions, and the other with 
the ratio of the sizes of small and large extra dimensions. In this paper, we discuss the impact 
of the radion fields on cosmology. We emphasize various aspects of radion physics such as radion 
coupling to the Standard Model fields, bare and dressed radion masses during inflation, dynamical 
stabilization of radions during and after inflation, radion decay life time and its late dominance 
in thermal history of the Universe as well as its quantum fluctuations during inflation. We argue 
that models where the radion plays the role of an inflaton or the inflaton is a brane scalar field, 
run into problems. We then present a successful inflation model with bulk scalar fields that seems 
to have all the desired properties. We also briefly discuss the possibility of radion as a cold dark 
matter candidate. 



I. INTRODUCTION 



The idea that the universe may have extra spatial dimensions has been under consider- 
ation for many years. The earliest thoughts date back to the works of Kaluza and Klein 
in the 1920 's. Such a possibility has been lately reinforced by the advent of string theory 
which is believed to be the most attractive candidate for a theory of gravity and gauge inter- 
actions lj, where the simplest realization requires at least 10 dimensions to be consistent. 
The extra dimensions are usually assumed to be compact with the sizes of the order of the 
fundamental Planck scale, l p = M~ l = (2.436 x 10 18 (GeV) -1 ~ 10~ 33 cm. However recent 
developments in this field have shown the prospects of having large compact dimen- 
sions 3|. The largest possible compact dimension could be as large as a few micrometers j^J. 
In the simplest incarnation it is believed that the Standard Model (SM) fields live in a four 
dimensional hypersurface or brane embedded in a higher dimensional space, which is known 
as the bulk, while gravity can propagate in all dimensions (for some earlier ideas see also 
Refs. p). 

In these models gravitational interactions appear to be stronger when probing distances 
below the size of the extra dimensions. The actual fundamental scale of gravity, M*, is no 
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longer the (now effective) 4 dimensional Planck scale, but it is rather related to this by 

M p = Ml +d -vo\- (1.1) 

where vol denotes the volume of the d compact spatial dimensions. Current measurements 
on the precision of the Newton's law at small distances impose an upper bound on the size 
of the extra dimensions at a level of 100 microns (|. For d — 2 it has been pointed out that 
the fundamental scale of quantum gravity could be as low as few O(TeV), which will have 
interesting implications for collider physics [zl^li]]. From this perspective these theories have 



an edge of being testable in future terrestrial and extra terrestrial experiments. The models 
with extra dimensions also provide an unique geometric insight to the naturalness problem 
in physics, which is one of the greatest reasons for studying various aspects of these theories 
in particle physics and in cosmology and astrophysics. Actually the stringent bounds on the 
size of the extra dimensions, or equivalently on the fundamental scale have been obtained 
from the astrophysical implications of the excited states of graviton, M* > 1000 TeV jlflj] . 
However very recent study reveals that such a bound can be weakened down to few TeV 
when the momenta associated with the extra dimensions are broken jll|. 



Another class of extra dimension models proposes that the SM degrees of freedom also 
g ate in some of the extra dimensions; examples of such studies can be found in Refs. [12I 
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. In this scenario, each SM particle is accompanied by a Kaluza 
Klein (KK) tower, with masses spaced by the inverse size of the extra dimensions in which 
they reside. Since there is no evidence of a KK tower associated with any SM particle 
yet, the size of such extra dimensions has to be smaller compared to the inverse of the 
electroweak scale, ~ TeV -1 . These models go by the name of universal extra dimension 
models (UED). Amongst the interesting features of these models is the possibility that 
they can be tested in near future in collider experiments through direct production of the 
excited KK states of the standard model particles Furthermore in six dimensions, 

such scenarios acquire the additional virtue that they can provide a natural explanation for 



the number of generations |l5l . Il6| , and an understanding of proton stability 
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despite the absence of any large 



with the lightness of the observed neutrino masses III 
mass scales in the theory. These models are therefore phenomenologically more sound than 
models where the standard model resides in 3+1 dimensional Minkowski brane. 

Cosmology of models with large extra dimensions, where only gravity pro pagate s in the 
bulk has received a great deal of attention, from building inflationary models 



25|,|26J to studying moduli problem 
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baryogenesis |29j. In this paper our main aim is to address similar cosmological issues in 
the class of models where both standard model particles and gravity propagate in the bulk 
(or the UED models). 

For simplicity we consider the n extra spatial dimensions where the standard model 
particles propagate, compactified on orbifold-ed circles with a common size r. The power 



law ruling of the gauge couplings in these m ode,s 113, Q and the need for keeping the 
theory from violating unitarity bounds |14lj . suggest that r should not be much larger than 
about 10-100 times the inverse fundamental scale (M" 1 ), i.e. M*r ~ O(10) - 0(100). 
Explaining the weakness of gravity in such theories (where M* of order of 10 TeV), requires 
that there must be large extra dimensions orthogonal to the 4+n already considered j^| (i.e. 
orthogonal to the SM brane). 

In the simplest scenario, we consider UED models comprising two sets of asymmetric 
extra spatial dimensions: n-small dimensions, of a size r just around a factor of ten to 
hundred times above the fundamental scale M" 1 , and p-larger spatial dimensions of size R. 
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Note that the SM fields are glued to 4 + n dimensions only, while gravity is propagating in 
the entire bulk. Therefore the effective four dimensional Planck mass M p is given by 

Ml = Ml +n+p r n R p . (1.2) 

For M* as low as few TeV, there should be at least two large extra spatial dimensions of sub 
millimeter size in order to yield the correct four dimensional Planck scale. 

In the above mentioned set-up, upon dimensional reduction, there are at least two radions. 
One describes the variations of the whole volume of the extra dimensions. We will call it the 
volume radion. The other, here called the shape mode, describes the relative expansion and 
contraction of the two sets of extra dimensions. In this paper we will take a phenomenological 
approach in determining the masses of the two radions, which could be either light or heavy 
compared to cosmological scale, which is usually taken to be the Hubble parameter. We 
will show that these radions play an extremely important role in determining the fate of the 
early Universe in UED models. 

The early universe cosmology is incomplete if it cannot reproduce the successes of the 
hot big bang model, e.g. synthesis of light elements at a scale ~ 1 MeV and their present 
day observed abundance. However this cannot be reached if the universe prior to this era 
were not already filled by the SM degrees of freedom, and provided there were no excess 
generation of entropy after big bang nucleosynthesis (BBN), which would otherwise dilute 
the initial abundance. This poses a challenge to any particle physics theory with large extra 
dimensions. In this regard UED models are no exception. They come with light scalar 
fields and therefore it is paramount that we study the dynamics and decay life time of 
these radions. The late oscillations and eventual decay of the radions can dominate thermal 
history of the universe, which can pose a serious threat to BBN. This is analogous to the 
string moduli problem. Therefore it is important that radions either decay before BBN, 
via their couplings to the SM fields, or they are sufficiently light to survive the age of the 
universe. If they survive late then they can be regarded as candidates for cold dark matter. 
In an ideal situation it would be possible to stabilize both the volume and the shape modes 
during inflation and/or right after the end of inflation. This is a major issue which we will 
deal in some detail in our paper. We argue that only the volume radion can be stabilized 
dynamically during inflation, because it always couple to the inflaton, whereas the shape 
mode remains uncoupled to it. Besides this, there is another important difference: the 
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volume radion governs the four dimensional Newton's constant, therefore it is important 
that we stabilize this radion before BBN. There are already stringent bounds on variations 
on the Newton's constant during BBN [3^ . On the other hand the shape mode can remain 
light and dynamically active field. We will find that it must also be stabilized before BBN. 
The shape mode can give rise to the running of the SM coupling constants. 

Furthermore, we also describe a successful inflationary scenario in UED model. We will 
discuss three possibilities, radion driven inflation, brane field driven inflation and the bulk 
field driven inflation and assess the relative merits of the different scenarios. 

Light scalars in cosmology are sometimes regarded as boon if they can sustain a flat 
potential with a large vacuum energy since they can then give rise to accelerated expand- 
ing phase of the universe, the inflation phase. They can also dilute all energy densities 
except the quantum fluctuations which can be imprinted on the cosmic microwave back- 
ground radiation. The radions in UED models could be responsible for generating adiabatic 



density fluctuations seen by COBE BOOMERANG |3j, MAXIMA |3j(, DASI 



candidate for inflaton, they can still be responsible for reheating the universe with the SM 
degrees of freedom and also converting its fluctuations into adiabatic modes. 

This paper is organized as follows. We begin by introducing our general set-up and 
performing a dimensional reduction of the gravitational part of the action in order to identify 
the radion fields. In the subsequent sections we devote our discussion on general aspects 
of the radion physics; we introduce the radion to matter couplings in section 3; we discuss 
initial conditions of the radion vacua before the onset of inflation in section 4. In section 
5, we comment on the problems of dealing with the radion (stabilization) potential and 
the radion mass, which seem to have an important effect on the dynamics of inflation and 
thermal history of the universe. Then we discuss the decay rate and life time of the radion 
fields. We begin section 6 with the analysis of the potential role of the radions in the 
early Universe. First we discuss the case when the radion acts as an inflaton. Then we 
extend our scenario by including a brane or a bulk inflaton. Then we discuss a possible role 
of radions to generate SM relativistic degrees of freedom along with the adiabatic density 
perturbations. We conclude our discussion by mentioning the possibility for the radions to 
become a candidate for the cold dark matter. 
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In this paper we propose that even if the radions are not the ideal 
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II. THE SET-UP, DIMENSIONAL REDUCTION AND THE TWO RADIONS 



Let us consider the following metric in 4 + d dimensions 

ds 2 = G ABdx A dx B = g iliV dx p 'dx v — hijdy l dy^ , (2-1) 

with A, B running over fi and i, where /i = 0, 1,2,3 labels our usual four dimensions and 
i,j — l,...,d are the total d extra dimensions. The coordinates y % take values in the interval 
[0, 1]. Note that since we have taken y % to be dimensionless, our hij has length dimension 
two and has L~ 2 . 

Also note that we are not considering the presence of vector-like connection A 1 ^ pieces, 
which are common in Kaluza Klein theories. This is because we are only interested in the 
zero mode part of the metric, after orbifolding the extra spatial dimensions. We have A^ odd 
under parity transformation and therefore it vanishes at the zero mode level. By reducing 
the 4 + d scalar curvature term, we obtain at zero mode level the effective action 

S = — fd^xd d yy/^GR[G] 

= _L J fa ^ XjjL [ R [ g ] _ \ d(i d » hi . _ -h^d,h i3 ■ h kl d»h k ^ , (2.2) 

where = M 2+d /2 defines the true fundamental gravity scale, while 1/k 2 = M 2 /2, 

defines the effective 4 dimensional Planck scale, which is related to the fundamental scale 
by a volume suppression, k 2 = k 2 ■ volo, see Eq. (jl.2j) . Here voIq is the stabilized volume of 
the extra space and corresponds to stabilized values for the h^. 

In order to obtain the 4 dimensional scalar curvature term in a canonical form, we have 
to perform a conformal transformation on the metric, 

designed to cancel the extra ^/h/vol coefficient of R[g] in Eq. ()2.2|) 

e 2v Vh/vok = 1 . 
The action in Eq. ()2.2|) is then transformed into 

S = 1 Jd 4 xV=g- {R[g} - ~<9^ d^hij + ^Wdphij ■ h kl d»h kl } . (2.5) 



(2.3) 



We take ip such that 



(2.4) 
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Next for the four dimensional part of the metric, g^, we will assume the standard Friedman- 
Robert son- Walker metric with a flat geometry i.e. 

= diag(l, -a(t), -a(t), -a(t)) , (2.6) 

for an isotropic and homogeneous (brane) universe, whereas we consider a diagonal form for 
the h part of the metric: 

n— times p— times 

h i:i «-> diag{h\, . . . , h\, h\, . . . , hi} , (2.7) 

where we have assumed n (p) of the g? extra dimensions to be small (large) with common 
physical sizes, hi(t) and /i 2 (^) respectively. Obviously d = n + p, and the physical volume 
of the extra space is given as 

volphys = Vh = h\h\ . 

At present, the physical sizes of the extra dimensions are given by hi t o = r and /i2,o — 
where the zero subscript stands for the stabilized values. Therefore the stabilized volume is 
given by vol® = r n R p , as stated in Eq. (jl.2j) . 

With the help of Eqs. (J2.7)) and (J2.5|) we obtain the non diagonal kinetic terms to be 

_L{ n(n + 2) (^) 2 + 2 „ p (^) +p(p + 2) (Ma) 2 } . ( , 8) 

We diagonalize this by introducing two physical radion fields 

"■(£)■ —^-(ss) ■ <M) 

This makes the kinetic terms in Eq. (|2.8|) diagonal and canonically normalized. Both the 
radions have a straightforward physical interpretation, a v is a sole function of \/7i, therefore 
it is related to the variation of the physical size of the volume, and often called a volume 
modulus (or volume radion). On the other hand <jj_ is related to the ratio of the sizes of 
both large and small dimensions, and we will often call it the shape mode. Furthermore, 
note that both the radion fields are set to be zero when stabilized (that is when hi and /12 
obtain the values r and R respectively). 

The effective four dimensional action for the gravity sector can then be written as 

S = Jd 4 x^-g {^- 2 R[g] + \{d,a v f + \{d^ i _f -U(a v ,a x )} . (2.10) 
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Note that in a particular case when n = 0, a± = and we are reduced to a single radion 
scenario of large extra dimensions. The last term in the above action represents the radion 
potential that must be introduced by hand to stabilize the extra dimensions. Its origin is 
largely unknown so far. In this paper we will only assume that there exists such a potential 
U(a v , <7_|_) and consider illustrative forms for it in our discussion. As we will discuss later on, 
the detailed form of U(a v ,a±) is important in understanding the dynamics of the radions 
during and after inflation. 



III. RADION COUPLINGS 



Given the above definitions for the radions, we now discuss how they couple to matter, 
both in the bulk and on the brane. We start by considering a bulk scalar field, <fi which 
latter on could be identified as the inflaton. The action, before performing the conformal 
transformation on the metric, goes as 



= J d 4 xd n+p y^/=g~Vh 



-G AB d A <pd B <p-v bulk {<p) 



(3.1) 

We now take only the zero mode level of the theory, therefore we write <p(x, y) = <p{x) / ^voIq, 
and properly scale the potential by introducing an effective four dimensional potential, 
V(<j>) = voIq ■ Vbuik((p)- By integrating out the extra dimensions and introducing the 
conformal transformation, we obtain an effective scalar action in the Einstein frame, 



Sd, = [ d 4 x\f^ge 2ip — - 

J VOlr 



(3.2) 



It is now straightforward to see that the last term gives rise to a coupling of <fi to a v due 
to the definition of ip given in Eq. ()2.4|) . This will also take care of the global exponential 
term, which cancels the y/h/volo piece of the action. Therefore we can simply write 

1 



Sc/, — J d x \f^g 
where the coupling constant is given by 



(3.3) 



a 



2(n + p) 



(3.4) 



' n + p + 2 

It is worth stressing that the a± radion does not couple to the bulk scalar field. Further- 
more, since the radion couplings appear due to the conformal transformation in Eq. ()2.3|) . 
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which depends only on a v , we can conclude the same is actually true for any bulk field. 
Therefore the a± radion will only couple to the brane matter, which we shall discuss now. 
Note that at this level we already notice an important difference in the dynamics of the two 
radion fields. 

Let us now consider a brane scalar field, x- The general action would be given as in 4 + n 
dimensions, where the SM fields propagate. This means that we have already integrated out 
extra p large dimensions which are harmless for the brane fields. Once more we stay at the 
zero mode level, we introduce a properly normalized zero mode by defining x — xl V^- We 
start with the effective action 



S x = I d A xy/=g [ -i 



n ^ 



r 



^d^X-Vix) 



(3.5) 



Notice that unlike the bulk field case, the induced metric does not contain the whole Vh, 
but only the piece that goes along the extra dimensions of the brane, h™. Therefore after 
conformally transforming the metric there will be a remnant of the h metric. Indeed by 
writing h™ = \/h K^ v and then, using the radion definitions in Eq. ()2.9|) in order to replace 
the /12 term, we obtain 

1 



S x = J d 4 x^ e ~^ /Mp e ~^ /Mp 
where the coupling constants are given by 



(3.6) 



^(n + p)(n + p + 2) 

7 = (3.8) 

We now see that the a±_ radion has a global coupling to both kinetic and potential terms of 
the brane fields, whereas extra global coupling of a v to both kinetic and potential terms also 
appear. Such couplings affect the definition of the canonical brane field, which are, however, 
nicely recovered around the minima of a Vj ±. 

As far as the other fields go, following a similar procedure, we find that the radion 
couplings to brane matter are given as: 

S gauge = ~ j rfxy/^jj e^W^ e ~l^/M p ^ ^ gW g vr ^ (3 9) 

for the zero mode gauge fields, and by 

= [ d 4 x^ e~ {a/2+ ^ /Mp e~^ /M - (-tyDtfail)) eZgV, (3.10) 



for the zero mode fermions. It is worth mentioning that if we take n = 0, as for the case 
of large extra dimensions, we find that a = (3, and thus the radion coupling to a massless 
gauge field vanishes Qj. This is however not the case when the matter fields propagate only 
in part of the extra space as shown in Eq. (J3.9j) . 



IV. INITIAL CONDITIONS BEFORE INFLATION 

The most natural scale for the universe to start in 4 + d dimensions is the fundamental 
scale. Therefore, by taking hi 2 (t = 0) ~ M" 1 as the initial values for the size of the extra 
dimensions, we obtain the initial vacuum expectation values (vevs) for a v and a± are of the 
order of the Planck scale: 



ln + p + 2, [M p \ fn ^ %f I np , (R 



Note that both the initial values for the radions are negative. This is due to the fact that 
the extra dimensions are relatively larger today (and so the volume). Physically, the extra 
dimensions have to grow, while both the radions approach to zero from negative values. 
For the sake of illustration we can take M* ~ O(TeV), and we find that indeed both the 
absolute values for the vevs, though n and p dependent, will have initial vevs quite large, 
e.g. a v (0) ~ O(10Af p ), a±{0) ~ O(10M p ). 

There is a couple of interesting points to be mentioned. First, the observed Universe is 
many orders of magnitude larger than M" 1 , and it is also certainly larger than even the 
largest p extra dimensions, which should be smaller than few micrometers in order to evade 
detection at small distance gravity experiments. Therefore an asymmetric expansion of the 
spatial dimensions is required. Besides the galaxy formation requires that there has to be 
perturbations seeded on large scales (larger than the size of the horizon). Inflation is the 
only paradigm which makes the universe large, flat, isotropic and homogeneous, if it began 
at M~ x . Even though it creates a homogeneous and isotropic universe, it cannot really 
beat the quantum fluctuations which are being stretched outside the horizon at the time of 
inflation. These quantum fluctuations when they reenter into our horizon late during the 
radiation dominated era provide the seed for the galaxy formation 37 1. 

However what remains to be seen is what could be the inflaton in our scenario. We 
should also remind the reader that at the moment fundamental theory provides no guidance 
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regarding the radion potential and, therefore the masses for a v , a± are not clearly fixed by 
our action Eq. (|2.10|) . This leads to a degree of speculation which we will discuss briefly. 

V. RADION MASS 

Although, from a naive point of view one would expect that the radion masses ought 
to be somehow related to the compactification scale, this seems to be less obvious when 
we recognize that physical radions in an asymmetric bulk are not solely connected to a 
particular extra dimension, but a combination of all. Moreover, since we have no knowledge 
about the origin and actual form of the stabilization potential, very little can be said about 
the radion masses without further assumptions. It has been conjectured that the radion 
mass can be as small as m r < M^/M p 38[ for M c = 1/r. This seems certainly true for the 
case where the stabilization potential comes from brane physics, however the mass can be 
much larger if bulk physics plays a role as we argue below. 

For the case where the stabilization potential is based on bulk physics there is a volume 
enhancement in the effective four dimensional theory. On pure dimensional grounds, we can 
see that the potential could be as large as M^M^, rather than Mf leading to much larger 
radion masses. For instance, in this case, we could write the radion potential roughly as 
U(a) < M^MpU(a/M p ), so that the upper bound on the radion mass comes naturally as 
m r ^ Another intuitive way of seeing this result is by noticing that masses for the 
radions jof , where % — v, _L, are invariant under dimensional reduction. The possibil- 
ity of having a large radion mass has also been pointed out in another example given in 
Ref. [33], where radion mass and compactification scale are exponentially related to each 
other, therefore reconciling the likelihood of large dimensions with large radion masses. 

As it has already been noticed 



38( a light mass radion can mediate long range interactions 
with a similar strength as gravity, therefore, the experiments testing the Newton's law at 
small distances p] impose a lower bound on such radion masses, which corresponds to the 
mass scale, m r > 10~ 3 eV. 
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A. A phenomenological guess on radion potential 

While we have no guidance from a fundamental theory regarding the nature of the radion 
potential, there are some cosmological constraints which could be used to guess its overall 
nature. First of all, we would like to have the minimum of the radion potential at a v = 
and cr_i_ = 0. Secondly, since the volume modulus, a v governs the magnitude of the Newton's 
constant, it must be stabilized before BBN, in order to be compatible with stringent bounds 
on the variation of the Newton's constant during BBN jscj. The shape mode must also be 
stabilized in order to evade constraints from long range interaction between the SM fields 
and late time variations of coupling constants. Since the stabilization crucially depends on 
the actual form of the stabilization potential, this provides some constraint on its form. 

The two extreme possibilities for the potentials are as follows: 

(i) The potential for the radions could be fairly flat, even for larger vevs, in which case a 
simple mass squared-like potential, e.g. U ~ \m 2 r o 2 r , where r = v, J_ can be used as 
a phenomenological model. As we will show this choice of potential can be extremely 
hard to stabilize. Besides the late oscillations of the radions will give rise to a moduli 
problem, for which there is no simple solution. 

(ii) This problem can be ameliorated provided that the moduli potential is steeper than 
mfof at large vevs, or in other words, the flatness of the potential is lifted. Let 
us consider for example a phenomenological potential that can give rise to a volume 
stabilization 

U(vol) ~ /i 2 M, 2+2(n+p) Mo - vol phys ) 2 , (5.1) 

where the overall scale has been taken in the general philosophy of using M* to adjust 
the dimensions of the potential. Here fi is assumed to be some mass scale naturally 
related to the volume, and so this is at most of order M^/Mp. By using Eq. (|1.2j) . we 
can rewrite the potential as 

U(vol) = m\Ml (1 - vol ph y S / W ) 2 = m 2 v M 2 (l - Vh/vol ) 2 , (5.2) 

where now the mass scale is m v = fiMp/M*, which is much larger than //, but yet 
smaller than M*. The above equation translates into a potential for the volume mod- 
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ulus with an exponential profile: 



U(a v ) ~ m 2 v M 2 p (l - e a ^ /Mp ) 2 . 



(5.3) 



Furthermore, note that the radion mass around the minimum becomes m v , therefore 
reinstating our arguments. Some interesting facts about the above potential are that 
(i) it is exponentially steep for a v > 0, which means that such a potential will not 
allow the volume to grow much larger than its stable size, (ii) On the other hand, 
the potential is fairly flat for large and negative vevs of a v , allowing the volume to 
easily expand. Let us stress here that we are simply using this potential for illustrative 
purposes. 

Needless to mention that the above arguments should also hold for the shape mode. 
However unlike the volume modulus it is extremely hard to make any speculation on the 
mass for the shape mode. The only obvious constraint comes from the long range interaction 
mediated by a±, which gives the lower bound m± > 1CT 3 eV, as already mentioned. In 
summary we can safely argue that the mass range for both the radions could be 



although there are models where the mass of the volume modulus radion could easily be 
close to the upper bound. 

VI. DECAY RATES OF THE TWO RADIONS 

Massive radions may decay into visible particles. Decay channels may vary depending 
on the actual mass. However, without loss of generality we can say that radions decay 
into photons, and possibly other SM particles. The couplings that are responsible for such 
processes are described in section three. In all cases, but when the SM fields are truly four 
dimensional (when n — 0), the decay would take place at tree level and, when evaluated 
around the stabilized radions they will always come with Planck suppressed couplings. The 
decay rate for any channel then goes as 



KT 3 eV < m r < M, 



(5.4) 




(6.1) 



192vrMj, ' 
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where the coupling constant a ~ 0(1) depends on the type of radion, specific decay channel 
and the number of dimensions, a is actually given in terms of the previously introduced 
coupling constants a, j3 and 7 (see section three). The life time of the radion is then given 
by 

r r « 1.2 x 10 7 yr (^-) . (6.2) 
V m r J 

It is illustrative to see that a radion of mass m r ~ 30 MeV has a life time r r « 10 11 yr, 
which is already about the age of the universe. Lighter radions would then be effectively 
stable, so they might be good candidates for the cold dark matter. For heavier masses 
the radions would decay into visible channels (photons and leptons) which might leave a 
detectable trace. For instance, for masses just around few MeV, radions can be produced 
not only in the early universe but also in supernova cores. Its later decay would inject hard 
photons of MeV energy at present time contributing to the diffuse gamma ray background 
flux, however, it is not too different from the one induced by KK graviton decays lOj, and 
so, no new bounds on the fundamental scale can be obtained from this effect. In contrast, 
a radion mass about a few TeV would decay into energetic SM particles just right at the 
BBN era (~ 1-10 sec), which may photo-disintegrate primordial nuclei. Obviously, a radion 
mass within MeVs up to just below TeVs should decay somewhere within BBN and present 
time. Radions with larger masses than 10 TeV, or so, would decay fast enough as to wash 
out any primordial radion abundance. They could potentially reheat the universe with the 
SM degrees of freedom. 

VII. RADION DRIVEN INFLATION 

The two radions evolve dynamically. It could be interesting to check whether by choosing 
an appropriate radion potential, one can drive inflation or not. Let us consider type I 
potentials for a v , a±. As already mentioned, irrespective of the masses m v , m±, the initial 
vevs are large in 4 dimensions, e.g. |o\,(0)| ~ | crj_ (0) | ~ O(10 M p ). At such large vevs 
the initial conditions are equivalent to chaotic initial conditions, and if there is no other 
source for cosmo logical constant, then inflation should take place along the line of double 

V ~\my v + \mlvl. (7.1) 
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Both the components give rise to the Hubble expansion, for which we obtain H pa 
10 max[m„,mi]. Using the mass window in Eq. (j5.4|) . we see that Hi n f could lie between 
O(10~ 2 eV) < Hi n f < 0(M*). Inflation will terminate when slow roll conditions are no 
longer valid, which happens in the chaotic model when a v , a± ~ y2M p . 

During inflation both the fields generate quantum fluctuations which leave the horizon. 
We can estimate the amplitude for their density perturbations as 

bp ^ max [m v , m±] M* 
p ~ M p < M P ' 

Note that the present constraint on the amplitude for the primordial density perturbations 
is given by one part in 10 5 , this is the COBE normalization verified by recent experiment 



such as WMAP |35l. l36l|. Unfortunately neither m v nor m± is sufficiently large to generate 
the right amplitude for the density perturbations in this radions driven inflation, as the 
bound on the RHS of the last equation indicates. 

An interesting possibility may arise if the total radion potential is given by a combination, 
such as 

V ~ m 2 v Ml (l - e <WM P ) 2 + \ m \ a \ , (7.3) 

where we assume a particular scenario with m v ~ 0(1 TeV), assuming for instance that 
there is some bulk physics which gives rise to the origin of volume stabilization, for e.g. flux 
quantization (4]]]. Let us also assume for the time being m± <C rn v . Note that here we get 
inflation from o v radion, because the potential has a flat direction for large and negative 
vevs, e.g. <J V (0) ~ — O(10M p ). Clearly for this case the exponential part on the potential 
is negligible, and the cosmological constant dominates, such that V ~ mr v M P . The Hubble 
expansion during inflation is given by Hi n f pa m v . Again it is easy to verify that the adiabatic 
density perturbations produced by a v during inflation will be way too small compared to 
10~ 5 , therefore one has to seek alternative ways of generating density perturbations. One 
such possibility is described below. 

Whenever there are more than one fields involved during inflation, the quantum fluctua- 
tions can be of two types, pure adiabatic fluctuations, which tells us about the perturbations 
along the final trajectory of the fields, and second one is the isocurvature fluctuations which 
designates the perturbations along the orthogonal direction of the trajectory. In general 
isocurvature fluctuations arise due to difference in pressure fluctuations and most impor- 
tantly isocurvature fluctuations feed the adiabatic fluctuations outside the horizon. 
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In our case a±_ being light compared to the Hubble expansion during inflation will also 
have homogeneous fluctuations. These fluctuations are actually isocurvature in nature. 

The amplitude for the isocurvature fluctuations generated by a± during inflation can be 
given by 42j 

V 1 ' 2 — — (7 4) 

In order to obtain the right amplitude for the density perturbations, a± ~ 10 5 x m v . If 
m v ~ 1 TeV, then the vev of a± ~ 10 8 GeV. However since the mass of a± is very light and 
the initial amplitude, i.e. |erj_| ~ O(10M p ) is quite large, it will be extremely hard for a± 
to roll down from M p ~ 10 18 GeV to the amplitude 10 8 GeV. Especially note that during 
inflation the radions are both slow rolling. Therefore the evolution of a± can be written as 

3Ha ± w -m\a ± . (7.5) 

Integrating this out we obtain 

(Xi(t)~ffi(0)e" (m ' w . (7.6) 

For m± <C H, it would take many e-foldings of inflation before the field could be settled down 
to its minimum. Anther way of saying this is that in order to obtain cr_|_ ~ 10 8 GeV from the 
initial value 10 M p , it would require between 76 < M e < 10 16 , for 1 TeV > m± > 10~ 3 eV, 
respectively, where Af e is the number of e-foldings. There is certainly a possibility that one 
of the radions could drive inflation while the other generates adequate density perturbations. 
If this can be realized then this will be the simplest model to be considered. 



A. Reheating & moduli problem 

Let us briefly discuss the reheating temperature due to the decay of radions. During 
coherent oscillations, the radions can decay into SM degrees of freedom. The couplings 
among radions and brane SM fields have been discussed in section 3. By looking at those 
expressions it is easy to see that the presence of a non trivial radion vacuum would have 
two effects: it can change the normalization of the wave functions through the global radion 
couplings to the kinetic terms; and it can affect the gauge and Yukawa couplings as well as 
the potential terms. Nevertheless the overall radion couplings remain Planck suppressed. 
Therefore for our naive estimations one may still use the results given in section 6. 
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Radion decay products must hadronize before BBN 43]. If however the radion decays 
much later than BBN then there will be two problems: (1) hadronizing the universe with 
SM relativistic degrees of freedom may not be complete, and (2) even if thermalization is 
completed well before BBN, the radion decay products with hard momentum can eventually 
destroy the light elements by reheating the universe. 

Using Eq. ()6.1|) the perturbative reheating temperature can be estimated by 

T r » 0.1 J M P T r = 1.28 x 1(T 7 GeV ■ — ^ . (7.7) 
v V Te V / 

Therefore, from our naive analysis, it seems that the only solution to this problem is that 
the radion o v must be heavier than about 500 TeV. Indeed for a mass at such value the 
reheating temperature is just about T r ~ 1 MeV. This also forbids the possibility of having 
a TeV fundamental scale in this scenario 1 . 

Let us discuss the associated moduli problem. In the above scenario a moduli problem 
will occur due to the coherent oscillations of a±. Note that a± oscillates with a frequency 
rri]_ when H ~ mj_. Depending on the initial amplitude of the coherent oscillations of <j± 
there can be adverse affects on cosmology. The energy density stored in p a± m m 2 ± a'j_ will 
redshift as in a matter dominated era and will eventually dominate the relativistic bath 
already created by the decay products of a v , provided p a « H 2 M 2 . However this requires 
the amplitude of <t_l oscillations to be as large as <J±(0) ~ M p . On the other hand, if the 
initial amplitude of the oscillations are sufficiently low, e.g. <7j_ < M p then it is certainly 
possible to avoid the radion dominance problem. This may occur if there is a large enough 
e-foldings of inflation. 



VIII. BRANE INFLATION AND ITS CONSEQUENCES 

In this section we consider a possibility where the inflaton potential arises from brane 
physics. The scalar field action is given by Eq. (|3.6jl . We suppose that a hypothetical scalar 
field x is living on the 4 + n dimensional brane and its potential is given by V(x)- The total 
action with the radions is given by 

1 Interesting possibility may arise if the reheating temperature is asserted to be the Hagedorn temperature, 
which comes out to be close to the mass of the radions, i.e. T r ~ 0.74 x m v < TeV, see [I^. 
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The natural scale for the effective V{x) in 4 dimensions is 

V( X ) < {M*r) n Mt ~ KTM* . (8.2) 

The numerical factor on the RHS arises when the compactification scale for n extra spatial 
dimensions is M c ~ lO^M*. Initial phase of inflation is governed by the exponential part 
exp [—((/? + a)a v + 7cr_|_)/M p ] K(x), since cr„(0), crj_(0) ~ —0(10 M p ), which exponentially 
enhances that contribution. 

In slow roll approximations, where we neglect the higher order time derivatives and we 
assume that the friction term proportional to the Hubble expansion dominates, the equations 
of motion for the fields and the Hubble parameter simplify a lot, and one gets 

3Ha v « P±^ e -W+^+'y^yMpV{x) , (8.3) 

IVlp 

3H X « - e -^y M "^- , (8.5) 
<9x 



2 



1 



3M p 2 



Let us assume that during inflation exp [— ((/3 + a)<r„ + 7cr_i_)/M p ] V"(x) ^> U(a v , a±). We 
can also assume that V(x) piece is dominated by a false vacuum and it is merely constant 
during inflation. However, note that a v and a± are present in the exponential part of the 
potential. The rolling of the radions will govern the initial phase of inflation in such a 
scenario. This is an example of assisted inflation, for which the scale factor follows 45] 

a{t) oc t K , k= 2 - . (8.7) 

[a + p) + Y 

Here we have assumed that V(x) is a flat potential, which means that x evolution is slower 
compared to a v , crj_. In our case the power law inflation lasts as long as \<r v \, \<yj_\ ~ M p . 
Once the vevs of both the radions become smaller than M p , then effective masses for the 
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radions become 

uu v p 

m 2 «^ + y»ZM (8 9) 

If V(x)/M% were greater than the bare masses for the radions, then an effective radion mass 
would become of order H in f ~ V(x)/M%. However, since V(x) ^ (M*r) n Mf, this means 
Hi n f is many orders of magnitude smaller than M*. For instance, for M* ~ 1 TeV and 
M*r ~ 10, we obtain an inflationary scale of order Hi n f < 10™ x 10~ 2 eV!. This also implies 
that the radions potential, U, should be extremely flat. Certainly H in f improves for larger 
vales of n and M*, but still it is hard to push the scale of inflation up to 0(TeV). Such low 
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scale of inflation has many problems to produce successful baryogenesis and BBN 

On the other hand if radion masses win over inflaton scale in Eqs. ()8.8)1 and (|8.9jl . then 
inflation ends right after a v , a± ~ M p . The radions start oscillations around their minimum 
and eventually decay into SM degrees of freedom. Note that the radion oscillations are large 
with the amplitudes of the order of M p . The late domination of radion can occur if there is 
a hierarchy in radion masses. The lighter one, if it does not decay completely, can cause a 
problem similar to the moduli problem as discussed before. 

IX. BULK INFLATION 

As a final possibility we investigate the bulk driven inflation. We assume that inflation 
arises from the bulk physics, e.g. a scalar field living in 4 + d dimensions, such that the 
action is given by 



S = J d'x^-g [~{d»a v y + -(<V ± ) 2 + -(«9,0) 2 - e-^'Vffl - U(a v , a ± ) 



(9.1) 



The four dimensional action is much simpler in this case. Note that after dimensional 
reduction the effective potential is V(4>) < M^M^, which can be relatively large compared 
to the former scenario where we assumed that inflaton potential as originating from the 
brane physics. Let us discuss this scenario in some detail. 

First of all note that the radion <j v only couples to the inflaton potential and not a±. 
Therefore a± does not get any effective mass correction through the inflaton potential, while 
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the volume modulus does. As we shall see this will play an interesting role in ameliorating 
the moduli, or late domination of radion in thermal history of the universe. 



A. Inflationary dynamics 

The equations of motion for the above action are given by 

a v + 3Ha v + ^^V^) - |^ = , (9.2) 
M p oa v 

dU 

a ± + 3H& ± - — = 0, (9.3) 
oa± 

dV 

(j) + 3H(f)-e- aav/Mp — = 0, (9.4) 
ocp 



1 



3M* 



2 2 2 y 1 



(9.5) 



In this particular scenario there are two phases of inflation possible. Initial phase of inflation 
is driven by the exponential part of the potential provided there is flat potential for V(4>). 
This phase is essentially driven by radion, and the subsequent phase of inflation is governed 
by V((f>) with a Hubble parameter H < M*. During inflation the slow roll equations are 
given to a good approximation by 

3Ha v » -^e a ^V(<P) + |^ , (9.6) 

Mp 0(J V 

dU 

SH& ± » — , (9.7) 
da v 

3H<j> w e -<*°v/M ? ( 9 8 ) 

! r e- aa " /Mp V(<f)) + U] . (9.9) 



H 2 



3MI 



The scale factor a(t) follows 

a(t) oc t (n+p+2)/n+p , (9.10) 

by assuming that V((f>) ^> C/ and evolution is extremely slow. Note that again we obtain 
a power law inflation. The power law inflation does not last for ever. Once a v ~ — M p , 
inflation is not supported along a v direction, rather inflation is supported by the inflaton 
potential V((f>). Also note that when the vev of a v becomes less than M p , we can expand 
the exponential term exp(— aa v /M p ). The effective mass for a v can be expressed as 

<e ff ^^ + ^^3aHl f . (9.11) 
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Where we have assumed H in f « V(cj))/M£, irrespective of the form of V(4>), provided of 
course that V((p) is sufficiently flat to support inflation. Further note that if the bare mass 
for a v is small compared to H in f ~ M*, then the Hubble induced mass term will dominate. 
Once Hubble induced mass term is switched on the radion o v will follow the evolution 

a v {t)^a vl e-^ff /ZH) \ (9.12) 

where a v , is the initial amplitude which is of order M p . The volume modulus therefore settles 
in its local minimum within one Hubble time. Note, however, that the global minimum for 
a v could be displaced during inflation. If such a displacement is large enough another phase 
of inflation could take place. That generally happens when the radion potential is quite 
flat 0]. We will discuss this issue in the next subsection. 

Further note that although the amplitude of a v can be reduced considerably during 
inflation, the same can not be true for the shape mode, a±. The other radion does not 
obtain a Hubble induced mass correction during inflation, because it does not couple to 
field at all at classical level. Therefore the amplitude of a± may not damp as much as we 
wish. Depending on its mass and its vev it might be possible that a± could be responsible 
for generating adiabatic density perturbations, which we will discuss later on. 

As noted earlier, we could easily get two subsequent bouts of inflation; initially driven 
by the radion a v and then by the inflaton field in the bulk. However number of e-foldings 
required for the structure formation depends on the exact thermal history of the universe. 

Especially 



2J that the 



Usually the number of e-foldings is much less than 60 if the scale of inflation is low^ 
if the scale of inflation is as low as H in f « M*, it was already shown in [ 
number of e-foldings required for structure formation would be 43, provided the universe 
reheats by T T ^ ~ 10 — 100 MeV. Therefore it is only the last 43 e-foldings of inflation which 
are important for the purpose of density perturbations. 

So far we have not specified the form of V(4>). However we have noticed that the upper 
bound on bulk inflation energy density is ~ M^M^. Therefore the second bout of inflation, 
which is mainly driven by the bulk field, 0, the scale of inflation, H in f ~ M*, which for 
M* ~ 0(TeV), is sufficiently low. It is easy to demonstrate that a single bulk inflaton 
model does not produce adequate density perturbations (see for instance the discussions 



in Ref. 



2l|), whereas a two field bulk inflation model, where inflation ends via a pha se 



transition similar to hybrid inflation scenario, can generate the COBE normalization 
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The effective four dimensional potential for two fields is given by 
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V(,,N) s (f )%X + £ ) V - »V + < > g) V + i™j#> , (9.13) 

where 0, iV are zero modes of the bulk fields. Note that the four dimensional couplings are 
Planck suppressed, which comes out naturally by integrating out the extra d dimensions. 
The naturalness condition requires the bulk couplings A, g ~ 0(1), and N Q ~ M*. Inflation 
in this scenario is supported by the false vacuum A 2 M^M^, because during inflation N settles 
in its local minimum, N = 0, and the rolling of field gives rise to the phase transition. 
However one has to ensure that the false vacuum contribution dominates over (l/2)m^4> 2 
term. In this respect it is plausible to give a small mass term for the field, <C M*. 
Note that in this scenario inflation ends when 

A fN M p 



(9.14) 

g V M* J y ' 

Note, if A ~ g, and, N ~ M#, we automatically get C ~ M„ . 

This potential was further studied numerically in Ref. 24|, and it was shown that in this 
scenario the phase transition which ends inflation happens extremely slowly. In the most 
natural case, Nq ~ M#, A, g ~ C(l), arid < 0.1M*, it is possible to have extremely large 
number of e-foldings of order, Af e ~ (9(10 6 ). 

The occurrence of large number of e-foldings can be understood due to the fact that 
the shape of the potential along the N direction is extremely flat due to extremely small 
coupling. The field iV rolls down extremely slowly and spends a considerable time near the 
top of the potential before rolling down to its global minimum N ~ v^M*, while the 4> field 
rolls slowly away from C towards its global minimum = 0. It can be shown that inflation 



ends in this picture when 
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2/(^)V(t)~mJ, (9.15) 

The phase transition could be made rapid, or the number of e-foldings during the phase 
transition can be made small, provided there is a hierarchy << N < M*. 



B. The Radion Problem 



Let us consider once more the potential of the radion-inflaton system, 

V(a v , 0) = e- a ^ M "V((P) + U(a v ) . 
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(9.16) 



Even though the global minimum do correspond to a v = 0, it is easy to see that during 
the phase of inflation driven by the inflaton potential the radion is trapped in a a false 
(effective) minimum. Indeed, if one takes the inflaton potential as effectively constant, the 
minimization condition for the total radion potential reads 

dU(<r v ) _ a ^ aav/Mp 



<9CX„ Mr 



Clearly, a v = is not longer a solution, unless V(<f>) were null. Thus, the minimum for 
the radion is displaced 46] to a position which strongly depends on the actual form of the 
radion potential. We should notice, however, that this effective minimum is always located 
on the positive range of o v values, which means that within a Hubble time the bulk has 
grown larger than its expected value (see the definition of the bulk radion in Eq. (|2.9|) ). 
This can be understood by noticing that the global minimum of a v is at zero, and the radion 
is originally coming down the potential from negative values, therefore, only for o v > 
one could have a positive slope for U(a v ), as it is needed to match the RHS of Eq. (j9.17j) . 
Moreover, this equation indicates that the effective minimum is located at a a v value where 
the slope of the radion potential catches up with a now suppressed inflaton potential. 

To clarify these points, let us take for instance the potential U(cr v ) = ^m 2 crl, where the 
radion mass is assumed to be lighter than the inflaton potential scale, Mi. Here Mi is taken 
to be the inflation scale as given by the inflaton potential alone (i.e. we use V(4>) = MpMf). 
In Figure 1 we have plotted the total radion potential for a particular case where aM; /m = 
10. There we have also depicted each term on the potential. Now one can easily sees why 
and where the minimum gets displaced. First of all the inflaton part of the potential drops 
very fast due to the exponential suppression of a growing radion, however, this is not fast 
enough as to catch with the radion potential which, being quite flat, grows much more 
slowly. Next point is that the exponential has no minimum, so the minimum appears only 
due to U (<7„), and therefore, such a minimum should be close to the a v value where both the 
potentials are almost equally relevant. Unfortunately, for the case we are considering, this 
happens when a v is about three times larger than M p . This effect is more dramatic when 
the hierarchy among m and Hi is larger. In such a case the inflaton potential in our picture 
is shifted to the right, which means that the minimum should be even at a larger o v value. 

This is troublesome for our previous picture of inflation: it means that, although the 
inflaton is slow rolling, its contribution to the Hubble expansion is being affected by the fast 
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FIG. 1: Total radion potential for aMj/m = 10 (continuous curve). Dashed curve corresponds to 
the radion mass potential whereas the dash-dotted curve plots the inflaton-radion coupling term. 
Notice the minimum of the total potential occurs close to the point where both terms are of same 
order. 

evolution of the exponential suppression that comes with a growing radion. Once both the 
terms in the potential are comparable we cannot neglect the radion contribution to inflation 
anymore, which means that a third phase of inflation would take place. In this new phase 
both fields contribute to expansion, however, now Hubble scale would be much smaller than 
Mj by some orders of magnitude. In fact, by using Eq. (|9.17|1 . we find that the equation for 
the minimum satisfies 

M P yY> 



Using t 
obtain 



lis expression one can evaluate the whole potential at the effective minimum and we 



V(a v ) = m 2 a v 



1 M P 
cr v + 



(9.18) 



2 a 

In other words, the scale of inflation should now be H in f ~ m. 

The system will remain in the false vacuum as long as the inflaton remains constant. 
Inflaton was, nevertheless, slow rolling right before the radion field settles down to the false 
vacuum. As inflation proceeds, the inflaton slowly takes smaller vevs, as does the energy 
density residing in its potential. For small values of the scale of the V(4>) potential, Mj, 
the minimum gets shifted to smaller values of a v . It would look like if the radion were 
slowly rolling down an effective potential given by Eq. ()9.18|) . By looking at the RHS of 
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FIG. 2: Same as before but now taking aMi/m = 0.1. Notice that here the minimum is consistent 
with <t v = 0. 

Eq. ()9.9|) we find that the exponential suppression of the inflaton potential does contribute 
for a very slow evolution of the inflaton field, which behaves as if it were effectively frozen. 
All these would result in a large period of inflation, that may last plenty of e-foldings. As 
inflation continues in 4D, the radion rolls down and the extra volume now shrinks. It ends 
only when the inflaton leaves the slow rolling phase, by then it is very likely that the radion 
would still have a large, almost Planckian vev. Therefore, just after inflation ends the radion 
would again start oscillating with an amplitude of order of the Planck mass. Therefore bulk 
inflation may not help to completely solve the radion problem. 

On the other hand, this problem does not appear if the radion mass were larger than 
the scale of the inflaton potential. Figure 2 plots the radion potential when aMj/m = 0.1. 
Now the whole picture changes dramatically, the radion mass dominates on the dynamics 
of expansion in a first phase. If the mass is large enough the radion quickly goes down the 
potential whereas the inflaton remains frozen. At some point, when the radion amplitude 
is quite small in the picture, the inflaton takes over and drives the radion towards the 
effective minimum, which is nicely just around a v = 0. From there on the evolution would 
be dominated by the inflaton potential and our discussion in the previous subsection applies. 
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FIG. 3: The total radion potential as seen when the radion stabilization potential U(a v ) becomes 
rather steep for positive values of a v . Here we have used aMj/m = 10. Notice that now the false 
minimum is closer to zero than before. 

C. A potential resolution to the radion problem 

Finally, let us consider what happens if the radion stabilization potential were similar to 
the one described in Eq. (|5.3|) : 

U(a v ) ~ m 2 v Ml (l - e <WA^ 2 _ (g 19 ) 

We use this potential and plot the total radion potential in Figure 3 for the same hierarchy 
of parameters used in Figure 1. Now, the stabilization potential grows exponentially for 
positive a v values and then catches up with the inflaton potential not far from a v = 0.5 Mp. 
In contrast to what happened in the case depicted in Figure 1, now the bulk volume is 
not allowed to grow too much beyond the stabilized value and this ameliorates the radion 
problem. Notice also that at the effective minimum the total potential is still dominated 
by the inflaton, therefore, it gives the features we were taking as granted in the previous 
subsections. Once the radion is settled at the effective minimum expansion would be driven 
only by the inflaton. This kind of potentials gives an idea of what one would desire for 
the actual stabilization potential. Therefore we conclude in a cheering note that indeed if 
the correct radion potential is found we would be able to stabilize the radion to its true 
minimum during inflation, and thereby solve the radion problem. 
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D. Reheating 



It is believed that the total energy density of the inflaton is transferred into radiation. 
The minimal requirement of cosmology is to have a thermal bath with a temperature more 
than 0(1) MeV in order to preserve the successes of BBN. The final reheat temperature 
depends on the decay rates T^n of the oscillating field 



Tr-O.l^ + r^Mp. (9.20) 

Note that in our model the bulk fields couple to effective 4 dimensions via Planck suppressed 
interactions with the SM fields on the observable brane. 

However the zero mode inflaton can reheat the brane and also the bulk. This possibility 
was already addressed in Ref. 22] for the case of gravitons. From Eq. (j9.13j) we found that 
both the fields obtain effective mass terms fn^ ~ fn^ oc \f2\N , if we assume A = g. Setting 
Nq ~ O(TeV), then both and N are kinematically allowed to decay into SM Higgs field, 



h. The decay rate can be estimated as follows |2 



/ 2 M 4 



<j>,N^hh ~ , 



32nM 2 m^ 



(9.21) 



where / is the coupling constant that we take of order one. If rh^ ~ rfiN ~ 0.1M*, then 
the reheat temperature is T r < 100 MeV, for M* ~ 100 TeV. Note that such a low reheat 
temperatures is a generic prediction if the inflaton field is living in the bulk. 

The next point is concerning the possible production of KK gravitons by the inflaton. The 
KK gravitons can be directly produced from the decay of KK modes of and N, through 
the channel n — > <ftiG n -i, where n, I are the KK numbers, and G is the KK graviton. On 
the other hand <p n modes can be produced via collision processes, such as 00 — > n 0_„, and 
similar reaction for N. The rate for exciting KK inflaton modes goes as 

Ml 

cr<t,4>-+cl> n <t>- n ~ ^TTi ' ( 9 - 22 ) 
p 

whereas the total decay rate of a KK scalar decaying via graviton emission is estimated as 

i 

12ttM2 



</) n ,total = 22 F<j>n-><l>eG n - e ~ 1 Ott- ]\ /f 2 ' (9.23) 



where = m 2 ^ + n 2 /R 2 . Note that the NN or 00 scattering rates for producing their KK 
counterparts are more Planck suppressed than the direct decay of 0, iV to the brane fields, 
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i.e. T^hh 3> <y<t><t>^4> n <l>- n x Ml ; where we have taken the number density of modes ~ Mf . 
So, we find that the zero modes of <fi and iV still prefer the Higgs as a final decay product. 
On the other hand, the excited KK inflaton or KK partners of other scalar fields present in 
the spectrum are extremely short lived. The heavier KK mode decays into the lighter KK 
modes plus gravitons, and eventually all the KK modes of the inflaton decay into the zero 
mode. In case where extra bulk fields are present the reheating of the bulk can be naturally 
avoided, if either those modes are as heavy as the inflaton, or the effective inflaton couplings 
to the bulk fields are smaller than the inflaton-brane interactions. 



X. DENSITY PERTURBATIONS 



A. Converting radion isocurvature fluctuations to adiabatic perturbations via ra- 
dion dominance 



We recognize that in a bulk inflation scenario the reheat temperature of the universe is 
small, T r ~ 1 — 10 MeV. The number of e-foldings required for the COBE normalization 
depends on thermal history of the universe after the end of inflation. For a low scale inflation 
and low reheating temperature the require number of e-foldings is small compared to the 



usual 60 e-foldings, it is roughly 43 e-foldings 



23, 



24j . As we discussed in our earlier sections, 



the last 43 e-foldings can easily arise during the slow phase transition for natural parameters, 
e.g. N ~ M*, ~ (10 _1 — 10~ 2 )M*. However as it was also pointed out in Ref. [24], the 
amplitude of the scalar density perturbations during the phase transition evolves gradually 
to become very small. Numerical studies suggest that for the above mentioned parameters 
the amplitude of the perturbations is way too small, e.g. V]^ < 10~ 15 [24]. The only 
way this situation could be ameliorated if there is a large hierarchy between the scales, e.g. 
Nq <C M*, rrifj, <C M*. The summary is that we do not gain on density perturbations by 
introducing two fields either. Nevertheless we should study other interesting options, such 
as radion induced adiabatic perturbations. 

Note that in our case a v can be dynamically stabilized during bulk inflation. However 
<j± cannot be dynamically driven to its minimum. As a consequence the actual dynamics of 
cr_i_ can be treated independently. If a± is sufficiently light compared to the Hubble expan- 
sion during inflation then quantum fluctuations along a± direction generates isocurvature 
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fluctuations, which could be transferred into adiabatic fluctuations when <tj_ decays while 
dominating the universe. However we have to keep in mind that a± must decay before 
BBN and it must dominate the energy density of the universe while decaying. In case a± 
dominates later on, then when it decays into SM degrees of freedom, it will inject significant 
entropy in the universe which will dilute the light nuclei created during BBN. Note that a± 
decays only via Planck suppressed interactions. This leaves us with a narrow window for 
the mass range of a±, which has to be large, e.g m a± = 500 TeV. 

Important issue to raise is what should be the amplitude of isocurvature density perturba- 
tions produced by a±. Let us for the time being concentrate on flat radion potential vr? a a\. 
As we discussed in section seven, see Eq. (17. 4|) . the amplitude of a± during that last 43 e- 
foldings of inflation must be 10 5 x M*. If we set an opportunistic value, m a± = 5 x 10 5 GeV 
and M* ~ 10 6 GeV, then <tj_ j j n / ~ 10 11 GeV, in order to satisfy the COBE normalization. 
Obviously, depending on the number of extra spatial dimensions, if the initial vev of a± at 
the onset of inflation is cr±(0) ~ — 10 M p , see section six, then following Eq. ()7.6|) we can 
evaluate the approximate number of e-foldings required to bring down its vev to a value 
when the the interesting scalar perturbation modes are leaving the horizon. The number of 
e-foldings comes out to be roughly M e ~ 3 x 10 5 . Note that here we are required to relax 
TeV as a fundamental scale. The main bound comes from the radion mass which has to 
decay before BBN in order to convert isocurvature perturbations to adiabatic ones. Also 
note that radion has to dominate the energy density before decaying. This we will discuss 
below. 

The conclusion is that if the fluctuations of a±_ generates V X J 2 ~ 10 -5 , then there has to 
be a large number of e-foldings of inflation. In this respect any chaotic inflation model can 
do the job, it could either be a single field inflation, e.g. V ~ m 2 cf) 2 with a large initial vev 
of (ft ~ 10 2 — 10 3 x M p , or two field hybrid inflation model with a very slow phase transition, 
see Eq. (j9~T3j) . 

The spectral index of the fluctuations can be determined by 

d\nV ai H m[ 2 m 2 
din A; Hf n{ 3# 2 nf 

If if i n f <^ H? ni , which obviously depends on the bulk inflaton potential, then there is a slight 
dependence on the spectral index given by n — 1 ~ (2/3)(m 2 L /if 2 lf ). On a generic ground one 
can say that if m cx < H ini ~ M*, then the spectral index is very close to Harrison Zeldovich 
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case, e.g. n ~ 1. This is one of the signatures of radion induced fluctuations. Note that the 
energy density of the radion is always subdominant during inflation, e.g. V > m 2 (7 a']_. 

However it is important that isocurvature perturbations stored in a± must be converted 
into the adiabatic ones. This can happen only if rn^ o\ dominates at later stages, but before 
BBN era. Therefore, the radion must be responsible for finally reheating the universe. 
However note that in our case this can happen soon after the end of inflation, because 
when m a± ~ Hit), the radion starts oscillating with an amplitude which is not far from 
10 8 GeV. Note that radion like the inflaton also couples to the SM degrees of freedom 
with a Planck suppressed interaction. As in our case both the inflaton and radion are 
oscillating simultaneously but with different frequencies, e.g. > m u± . However inflaton 
being heavier will preferably decay before the radion and towards the last stage, if the 
radion oscillation could dominate the inflaton decay products, then, it is possible to convert 
isocurvature fluctuations into the adiabatic ones. Otherwise there will be a contribution 
from isocurvature fluctuations also. The total curvature perturbations is given by the sum 
of 

V 1 * 2 ~ fpg + (1 - f)pW , (10.2) 
where / is the fraction of energy density stored in the radion, it is given by 

^3?fev (1 ° 3) 
where f2 7 is the energy density stored in the decay products of the inflaton. In our case 
this ratio could be small if f2 7 > fi CT± , especially since VJnf m^ ± aj_ and both the fields 
have a similar decay rate to the SM fields. However if we imagine that not all the inflaton 
energy density goes into the brane degrees of freedom, suppose the inflaton couples strongly 
to the bulk degrees of freedom, then it is quite possible that £l± 3> f2 7 . Note that even 
though there is a larger relativistic brane degrees of freedom, but the bulk volume is also 
large. Therefore it is possible to imagine that a branching ratio between brane and bulk is 
suppressed by the bulk volume factor. 

The conclusion is that it is indeed possible to generate adiabatic density perturbations 
from the isocurvature fluctuations of a±. However a realistic scenario requires larger funda- 
mental scale than TeV, e.g. M* > 500 TeV. 
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B. Converting radion isocurvature fluctuations to adiabatic perturbations via fluc- 
tuating inflaton coupling 



There is an alternative mechanism for converting isocurvature fluctuations of a± to adia- 
batic perturbations from fluctuating inflaton coupling. Note that in our case a v has a natural 
coupling to the inflaton while a± is an independent field which only couples SM degrees of 
freedom. Suppose the a v has not settled down to its global minimum during inflation, which 
may occur if a v has a flatter potential for both positive and negative vevs. In such a case 
the inflaton coupling to the SM fields will fluctuate when the inflaton decays due to the 
presence of the term e~ acTv ^ Mp V in the action. During inflation a v has fluctuations, and 
if it is lighter compared to H in f then it also has a non-zero vev. Further the above term 
can give rise to fluctuating inflaton coupling or fluctuating inflaton mass. In either case 
the perturbations can be transferred to the SM degrees of freedom during the decay of the 
inflaton. This can be understood in two steps; first note that fluctuating coupling or mass 
gives rise to fluctuations in the decay rate of the inflaton and second, fluctuating decay rate 
gives rise to fluctuations in the reheat temperature of the universe. The fluctuations in T 
can be translated into fluctuations in the energy density of a thermal bath with 47 1 

= (10-4) 

The factor 2/3 appears due to red-shift of the modes during the decay of the inflaton whose 
energy still dominates. It has already been proven that the curvature perturbation can be 

1 /2 

fed by the fluctuations in the decay rate of the inflaton V/ = —(1/6)ST/T 

that in this case, the tilt in the spectral index is given by Eq. (jlO.lj) . 

In our case, fluctuations in the decay rate is given by 

ST 5a v ^ H ini 
r a v 2na v 

In order to match the observed COBE normalization, the vev of a v during inflation must 
satisfy a v ~ 10 5 ifj n f ~ lC^M*. Since we know the initial vev of o v at the onset of inflation 
has to be around cx„(0) ~ 10 M p , therefore one requires at least Af e ~ O(10 4 — 10 5 ) number 
of e-foldings to bring down the vev of a v to appreciable value in order to generate the right 
amplitude for the density perturbations. 

In reality both the radions can play roles simultaneously; fluctuations in a v can give rise 
to fluctuations in reheat temperature of the universe and late decay of cr_i_ can transfer its 
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isocurvature perturbations to adiabatic ones. However the dominant one will prevail and it 
seems that both requires large number of e-foldings of inflation, which can be obtained for 
the bulk inflation only, see Eq. (|9.13|) . 



XI. RADION AS A CANDIDATE FOR COLD DARK MATTER 



Since the class of models we are discussing do not have the conventional dark matter 
candidates such as the axion or the neutralino, we like to explore other particles as CDM. 
As we have noticed a successful cosmological scenario seems to point to a reheat temperature 
of the universe which is much smaller than the cut-off scale M*. Therefore in order to have 
a stable KK dark matter 49J, it is necessary to excite them non-thermally. Nevertheless we 



should also keep in mind that we do not understand thermalization of the inflaton decay 
very well and it could be that for some reason the reheat temperature is close to T r < M* 
in which case KK modes would have no trouble getting excited and play the role of dark 
matter. 

There however exists another interesting possibility that the radion oscillations can also 
act as a cold dark matter provided the radions are absolutely stable and light. According 
to Eq. ()6.2|) this happens if the radion mass is lighter than 30 MeV. Thus, if we believe 
in a radion potential that traps the radion in the minimum during inflation, then once 
m Uv ~ H{t) <z 30MeV or m a± ~ H(t) <, 30MeV, the radions start oscillating respectively 
in their minimum with an amplitude of order H(t). These are coherent oscillations which 
give rise to an effective pressureless fluid and the density stored in their coherent oscillations 
can give rise to a cold dark matter. This issue will be dealt separately in near future. 



XII. CONCLUSION 



In this paper we have studied some aspects of the cosmology of models with universal 
extra dimensions, where SM particles live on a 4 + n dimensional brane, whereas gravity 
resides in 4 + n + p dimensions (with p dimensions being much larger than n dimensions). 
These models generically have two radion fields in the effective four dimensional theory. The 
two radions can be identified as (i) The volume radion, associated to the variation of the 
whole volume of the extra space; and (ii) the shape mode, produced by the variations in the 
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hierarchy among the n small and the p large extra dimensions in the theory. 

Our study shows the difficulties that the radions introduce in the description of early 
universe cosmology. In exploring this we consider two extreme forms of the radion potential. 
If the radion stabilization potential is rather flat, almost all scenarios of inflation suffer from 
a radion problem analogous to the well known moduli problem. We note that, although 
the radions themselves can drive inflation, the picture seems hard pressed to reconcile with 
the correct amplitude for the the primordial density perturbations and, typically inflation 
ends leaving the radions to oscillate with large amplitudes. We then focus on brane induced 
inflation and note that inflation occurs at relatively low scale and the inflaton energy density 
is typically too small (about 0(eV) A ) to drive a successful inflation. On the other hand, 
a bulk inflaton can help to settle down the volume radion to an effective minimum within 
a Hubble time and seems more promising from the point of view of density perturbations. 
Nevertheless, the effective minimum for the radion potential is usually located at large 
values of the radion field, which means that the bulk inflaton makes the volume of the extra 
dimensions to grow much beyond the expected stabilized size. If inflation is supported by 
both the volume radion and the inflaton, while the inflaton keeps rolling down the volume 
shrinks. In contrast, the shape mode gets decoupled and therefore remains frozen during 
the entire inflationary phase. By the end of inflation it is very likely that both the radion 
modes will end with large amplitudes, unless the universe passes through an extremely large 
number of e-foldings. 

Most of the above mentioned problems get resolved if the radion mass is assumed to be 
large. This allows the radion to decay fast into SM particles while oscillating at late stages. 
With a bulk inflaton, the radion problem can disappear completely for the volume radion 
which is now rapidly driven to its actual minimum and trapped there by the dynamics of 
inflation. For this it is mandatory that the flatness of the volume modulus is lifted by some 
brane or bulk physics. We have explored phenomenologically this possibility. 

Finally we concentrated on the radion induced isocurvature fluctuations. The observa- 
tions demand that we convert this isocurvature fluctuations to the adiabatic ones. It is quite 
possible in our case because the radions can decay into the SM degrees of freedom. However 
this again would require the radion masses to be fairly large > 10 TeV at least, so they 
decay before BBN. Another interesting possibility which we have noticed is that the inflaton 
coupling to the SM matter can fluctuate, because of the couplings being governed by the 
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vevs of the radions. These fluctuations induce fluctuation in the coupling and therefore in 
the reheat temperature of the Universe. This could be an interesting avenue to obtain the 
desired amplitude for the density perturbations. This particular scenario does not demand 
large masses for the radions, but it certainly requires fairly intermediate scale vevs for the ra- 
dion mass, e.g. ~ 10 5 M*. This can be obtained easily if there is large e- foldings of inflation, 
which occurs very naturally in the bulk driven inflation scenario. 

The punch line of this paper can be summarized as follows; for (9(100) TeV scale gravity 
the bulk driven inflation along with lifted potential for the volume radion can give rise to a 
successful cosmology. 
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